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The energy loss of a fast parton scattering elastically in a weakly cou- 
pled quark-gluon plasma is formulated as an initial value problem. The 
approach is designed to study an unstable plasma, but it reproduces the 
well known result in the case of an equilibrium plasma. Contributions to 
the energy loss due to unstable modes are shown to exponentially grow in 
time. An unstable two-stream system is considered as an example. 

1. Introduction 

When a highly energetic parton travels through the quark-gluon plasma 
(QGP), it losses its energy due to elastic interactions with plasma con- 
stituents. This is the so-called collisional energy loss which for the equi- 
librium QGP is well understood, see the review [I] and the handbook [2]. 
The quark-gluon plasma produced in relativistic heavy-ion collisions, how- 
ever, reaches the state of local equilibrium only after a short but finite time 

* Presented by K. Deja at the HIC-for-FAIR Workshop & XXVIII Max-Born Sympo- 
sium 'Three Days on Quarkyonic Island', Wroclaw, Poland, May 18-21, 2011. 



(1) 



2 



E-Loss- Wroclaw printed on November 22, 2011 



interval, and during this period the momentum distribution of plasma par- 
tons is anisotropic. Collisional energy loss has been computed for such an 
anisotropic QGP [3]. However, the plasma with anisotropic momentum dis- 
tribution is unstable due to the chromomagnetic modes (for a review see 
[1]), and the fact that unstable systems are explicitly time dependent as 
unstable modes exponentially grow in time has not been taken into account 
in the study [3| . 

Our aim is to formulate an approach where the energy loss is found as the 
solution of an initial value problem. The parton is treated as an energetic 
classical particle with SU(3) color charge. For the equilibrium plasma we 
recover the known results but for the unstable plasma the energy loss is 
shown to strongly depend on time. Our approach to the energy-loss problem 
is similar to the method used earlier to study the momentum broadening q of 
a fast parton in anisotropic plasma [5]. Analogous methods have also been 
used to study the spectra of chromodynamic fluctuations of an unstable 
plasma [B] which, in particular, are responsible for the collisions integrals of 
transport equations [7j. 

Throughout the paper we use the natural system of units with c = h = 
kB = ^ and the signature of our metric tensor is — , — , — ). 

2. General formula 

We consider a classical parton which moves across a quark-gluon plasma. 
Its motion is described by the Wong equations [8] 

u^{r), (1) 
gQ'^{r)Fr{x{T))u,{T), (2) 
-5r%(TX(x(r))Qe(r), (3) 

where r, x^(r), u^{t) and p^{t) are, respectively, the parton's proper time, 
its trajectory, four-velocity and four-momentum; Fa'^ and Aa denote the 
chromodynamic field strength tensor and four-potential along the parton's 
trajectory and is the classical color charge of the parton; g is the coupling 
constant and Og = q'^/'^tt is assumed to be small. We also assume that the 
potential vanishes along the parton's trajectory i.e. our gauge condition 
is Uij_{t) Aa[x{T)) = 0. Consequently, due to Eq. ([s]) the classical parton's 
charge Qdj) is constant within the chosen gauge. 

The energy loss is given directly by Eq. ^ with = 0. Using the time 
t = 'jT instead of the proper time r and replacing the strength tensor Fa'^ 



dr 

dr 
dQajr) 
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by the chromoelectric Ea(f,r) and chromomagnetic Ba(t,r) fields, Eq. ^ 
gives 

^=5g'^E,(t,r(t)).v, (4) 

where v is the parton's velocity. Since we consider a parton which is very 
energetic, v is assumed to be constant and = 1. Introducing the current 
generated by the parton r) = gQ'^^r5'^^\Y — vt), Eq. M| gives 



dE{t) 
dt 



J d^rEa{t,r)-Ut,r). (5) 



Since we deal with an initial value problem, we apply to the field and 
current not the usual Fourier transformation but the so-called one-sided 
Fourier transformation defined as 

poo p 

/(w,k) = y dty (i3re*(^*-'^-^)/(i,r), (6) 

with the inverse 

where the real parameter cr > is chosen is such a way that the integral 
over uj is taken along a straight line in the complex cj— plane, parallel to the 
real axis, above all singularities of /(a;,k). Using Eqs. ^ and Q, Eq. ^ 
can be rewritten: 

dE{t) ^„ /■°°+^'^ duj f d^k 



dt 



where a) = k • v. 

The next step is to compute the chromoelectric field E^. Applying the 
one-sided Fourier transformation to the linearized Yang-Mills equations we 
get the set of equations familiar from electrodynamics 

ik'e'^{uj,k)Ei{u},'k) = /)a(w,k), ik'Bi{u,k) = , 

ie'^'k^E'^iu, k) = iojBiioj, k) + SS,(k), (9) 
ie*^-^fc^i?„^(t^,k) = j:iu,k)-iu;e'^ico,k)Ei{io,k)-DU^), 

where pa is the color-charge density, the fields with the index are ini- 
tial values; the chromoelectric induction is expressed as D^(a;,k) = 
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e*-'(a;, k)£'a(w, k) with e^^{uj,\i) being chromodielectric tensor which carries 
all information about the medium. For an anisotropic plasma it equals 



(r_ 

2oj 



d?p df{p) 
(27r)3 a; - k • V + iO+ dp'' 



where /(p) is the momentum distribution of plasma constituents. The color 
indices a, h are dropped as e(a;, k) is a unit matrix in color space. 

Although equations ^ strongly resemble those of electrodynamics, the 
gluon contribution to the color charge density pa and color current ja, which 
is a genuine nonAbelian effect, is fully incorporated in these equations. 

Using Eq. S, the field £^*(a;,k) is found to be 



where 



i(S-i)^^(.;,k) u;j^(a.,k) + e^'''k^B\,{\^) - ujB{^[V) , (10) 



S*-'(a;, k) = -k^5*J + + u?s^'\u}, k). 



:ii) 



Substituting the expression (10) into Eq. ([8j), we get the final formula 



dE{t) 
dt 



oo+ia 



d^k 



2m J (27r)' 



'^(^-'^)*(S-i)^^'(w,k) (12) 



ZLjgQ'^v^ 



UJ — UJ 



+ e^-^'fc^<(k)-c^I?^o,(k) 



The integral over oj is controlled by the poles of the matrix S~^(a;, k) which 
determine the collective modes in the system. Equivalently, these modes are 
found as solutions of the equation det[S(a;, k)] = 0. Equation (12) needs to 
be treated differently for stable and for unstable systems. 



3. Stable systems 

When the plasma is stable, all modes are damped and the poles of 
S~^(cj,k) are located in the lower half-plane of complex w. Consequently, 
the contributions to the energy loss corresponding to the poles of S^^(a;, k) 
exponentially decay in time. The only stationary contribution is given by 
the pole a; = a) = k • V which comes from the current ja(<^, k). Therefore, 



the terms in Eq. (12), which depend on the initial values of the fields, can 



be neglected, and the energy loss of a fast parton in a stable plasma is 

^ = -^g'Cnvv| ^ O (S-i)^^(^, k), (13) 
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where the bar indicates that averaging over parton's color state has been 
performed, and the factor Cr is 4/3 for a quark and 3 for a gluon. 

When the plasma is isotropic, the dielectric tensor can be expressed in a 
standard way through its longitudinal[eL(a;, k)) and transverse (^^('^jk)) 
components and the matrix S*-' (w, k) ( 11 ) can be easily inverted. The energy 



loss (13) then equals 



dE 
~dt 



d^k 



(2^)3 k2 



2, ,2 



+ 



kV 



■ UJ 



(14) 



which corresponds to the standard energy loss due to soft collisions [2]. 



4. Unstable systems 

When the plasma is unstable, the matrix S^^(a;,k) contains poles in 
the upper half-plane of complex u, and the contributions to the energy loss 
from these poles grow exponentially in time. After a sufficiently long time 
the parton's energy loss will be dominated by the fastest unstable mode. 



For an unstable plasma, the terms in Eq. (12), which depend on the initial 



values of fields D and B, cannot be neglected, as these terms are amplified 
by a factor exponentially growing in time. Using Eq. ^ the initial values 
can be computed as 



^Oa(k) 



^Oa(k) 

duj 1 
27r Lu 



k X Ea(a;,k))' = -igQ^e'^^k^ {T.-^f\u),\^)i 



(15) 
(16) 



ton's color as before we obtain 

= q'Crv'v' / — 

at .l-oo+ia 



Substituting Eqs. (|15|) and (|16|) into Eq. (12) and averaging over the par- 

d^k 



(27r)S 



<'^-'^^\j:-y^{uj,k) (17) 



n 



UJ 



■ UJ 



Equation ( 17 ) gives the energy loss of a fast parton flying across the unstable 



plasma. The key point is that the presence of an unstable mode will produce 
a time dependent exponential growth of the form e^™'^^'')*. 



5. Two-stream system 

In order to calculate the energy loss, one must invert the matrix S*-' (a;, k) 
defined by Eq. ( 11 ) to substitute the resulting expression into Eq. ([TtI). For 
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a general anisotropic system this is a tedious calculation, and therefore we 
consider for simplicity the example of the two-stream system which has un- 
stable longitudinal electric modes. We assume that the chromodynamic field 
is dominated by the longitudinal chromoelectric field and take B(w, k) = 
and E(a;,k) = k(k • E(i:i;, k))/k^. Then, S*-^(a;,k) is trivially inverted as 



1 V 



(18) 



and Eq. (|17|) simplifies to 

-oo+ia^J (27r)3w2g^(a;,k)k2 



dE{t) 



dt 



+ - 

LO — OJ OJ 



(19) 



Eq. (19) gives a non-zero energy loss in the vacuum limit when sl — )• 1- 
Therefore, we subtract from the formula (19) the vacuum contribution, or 
equivalently we replace I/el by I/el — 1- 

The next step is to calculate £^(0;, k). With the distribution function of 
the two-stream system in the form 



/(p) = (27r)3n[5(3)(p-q) + <5{3)(p + q^ 

where n is the effective parton density in a single stream, one finds [6] 

_(k)) {u + w+(k)) (w - w_(k)) (w + a;_(k)) 



eL(tj,k) 



[UJ — OJ^ 



(20) 



(21) 



where u = q/i?q is the stream velocity, fi^ = g^n/2Ec^ is a parameter 
analogous to the Debye mass squared, and iba;-|-(k) are the four roots of the 
dispersion equation ei(a;,k) = which are 



4(k) 



k2 



k2(k-u)2+/i2(k2_(k.u)2) 



(22) 



It is easy to see that < i^+(k) G M for any k. For k2(k • u)^ > 2^^{k^ - 
(k • u)2), the minus mode is also stable, < i^-(k) G M, but for k • u 7^ 
and k^ (k • u)2 < 2/2^ (k^ — (k • u)^) one finds that uj- (k) is imaginary which 
is the well-known two-stream electric instability. 

Strictly speaking, the stream velocity u given by the distribution func- 
tion (20) equals the speed of light. However, the distribution (20) should 
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Fig. 1. The energy loss per unit length as a function of time for a parton flying 
along the direction of streams. The lowest (green) curve corresponds to fcmax — 20, 
the middle (red) one to /cmax = 50 and the most upper (blue) curve to /cmax = 100. 



be treated as an idealization of a two-bump distribution with bumps of fi- 
nite width. Then, the momenta of all partons are not exactly parallel or 
antiparallel and the velocity u, which enters Eqs. (21, 22), obeys < 1. 



Equations ( 19 ) and ( 21 ) determine the energy loss of a parton in the two- 



stream system. The integral over cj can be computed analytically as it deter- 
mined by the six poles of the integrand located at cj = ±a;_|_(k), iba;_(k), w 
and 0. The remaining integral over k must be done numerically. 

Performing the calculations we have redefined all dimensional quantities 
by multiplying them by the appropriate power of ^ to obtain dimensionless 
variables. We have chosen the following values of the parameters: g = 
1, |v| = 1, |u| = 0.9, Cr = 3. To take the integral over k, cylindrical 
coordinates with the axis z along the stream velocity u have been used. 
Since the integral appears to be divergent, it has been taken over a finite 
domain such that — /cmax ^ < femax and < fcy ^ ^max- In Fig. [l] we 
show the parton's energy loss per unit length as a function of time for a 
parton flying along the streams for fcmax = 20, 50, 100. The energy loss 
strongly oscillates with amplitude growing in time. However, the results 
depend on the cut-off parameter /cmax: the magnitude of energy loss grows 
and the period of oscillations shrinks when /cmax increases. This ultraviolet 
sensitivity of our results is not very surprising, as our approach is fully 
classical. In the case of equilibrium (stable) plasma, the energy loss due 
to soft interactions diverges logarithmically with /cmax [I]- The divergence 
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signals a necessity to combine the classical contribution to the energy loss 
at small wave vectors with the quantum contribution at higher ones. A 
quantum approach to the parton energy loss in unstable plasma needs to 
be developed. 

Although the two-stream system is not directly relevant to the quark- 
gluon plasma which is produced in relativistic heavy-ion collisions, let us 
comment on the potential interest of our study. Choosing /i = 200 MeV, 
which is a rough estimate of Debye mass in the QGP from nuclear collisions, 
the units on the horizontal and vertical axis in Fig. [T] are, respectively, fm 
and 200 MeV/fm. Since the parton's energy loss in the QGP observed in 
relativistic heavy-ion collisions is of order 1 GeV/fm, the magnitude of the 
energy loss discussed here is certainly of phenomenological interest. 

6. Conclusions 

We have developed a formalism where the energy loss of a fast parton 
in a plasma medium is found as the solution of initial value problem. The 
formalism allows one to obtain the energy loss in unstable plasma where 
some modes exponentially grow in time. In the case of stable plasma, one 
reproduces correctly the standard energy-loss formula. As an example of 
an unstable system we have studied a two-stream system. The energy loss 
per unit length is not constant, as in an equilibrium plasma, but it exhibits 
strong time dependence. 
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